3D Anderson transition for two electrons in 2D 
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It is shown that the Coulomb interaction can lead to delocalization of two electron states in two- 
dimensional (2D) disordered potential in a way similar to the Anderson transition in three dimensions 
(3D). At fixed disorder strength the localized phase corresponds to low electron density and large 
value of parameter Vs . 
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Contrary to the well established theoretical result pi , 
according to which noninteracting electrons are always 
localized in 2D disordered potential, the pioneering ex- 
periment by Kravchenko et al. |^ demonstrated the ex- 
istence of metal-insulator transition for real interacting 
electrons in 2D. The ensemble of experimental data ob- 
tained by different groups [^^ clearly indicates the im- 
portant role played by interaction. In the majority of 
experiments the Coulomb energy of electron - electron 
interaction E^,, is significantly larger than the Fermi en- 
ergy Ep , estimated for noninteracting electron gas in ab- 
sence of disorder. The ratio of these energies is character- 
ized by the dimensionless parameter r^ = l/y^7m7a|j ~ 
Eee/Ep, where Ug is the electron density in 2D, and 
a*g — h eo/m*e'^, m* , eg are the effective Bohr radius, 
electron mass and dielectric constant respectively. Such 
large r^ values as 10 - 30 have been reached experimen- 
tally M-M ■ At these rg the electrons are located far from 
each other and it is natural to assume that in this regime 
the interaction effects will be dominated by pair interac- 
tion. The important role of the residual two-body inter- 
action is also clear from the fact that in the Hartree-Fock 
(mean field) approximation the problem is again reduced 
to the one-particle 2D disordered potential with localized 
eigenstates |l[. 

The problem of two electrons interacting in the local- 
ized phase is rather nontrivial. Indeed, recently it has 
been shown that a short range repulsive/attractive in- 
teraction between two particles can destroy one-particle 
localization and lead to creation of pairs propagating on a 
distance much larger than their size p-[l3| . The pair size 
is of the order of one-particle localization length Zi . Inside 
this length the collisions between particles destroy the 
quantum interference that results in their coherent prop- 
agation on a distance Ic ^ h- The important point is 
that only pairs can propagate on a large distance. Indeed, 
the particles separated by a distance R^ li have expo- 
nentially small overlap, the interaction between them is 
weak and such states are localized as in the noninteract- 
ing case. According to the theoretical estimates |g,^,|3| 
in 2D the localization length Ic grows exponentially with 
li according to the relation \n{lc/li) ~ k > 1. Here 
K ^ ^2P2, where r2 ~ U'^/{Vlf) is the interaction in- 



duced transition rate between localized states in e.g. 2D 
Anderson model, p2 ~ ^i/^ is the density of two-particle 
states directly coupled by interaction, V is the hopping 
between nearest sites, U is on (nearest) site interaction, 
and energy is taken in the middle of the band. In a sense 
the above estimate is similar to the case of one-particle 
localization in 2D where In^i ^ kp£ ^ (V/W)'^ and the 
product of the Fermi wave vector kp on mean free path 
£ is proportional to a local diffusion rate [Q; W is the 
strength of on site disorder. Indeed, in the same manner 
the interaction induced diffusion rate of a pair is given by 
Z?2 ~ I1T2 ~ k/II cx In/c- According to the above esti- 
mates Ic should vary smoothly with the effective interac- 
tion strength characterized by the dimensionless parame- 
ter K. However, this consideration is valid only for a short 
range interaction while the analysis of the long range 
Coulomb interaction requires a separate study. The in- 
vestigation of this case is also dictated by the experiments 
|p|-p| where the electrons are not screened and are located 
far from each other (r^ 3> 1). On a qualitative grounds 
one can expect that the effect of Coulomb interaction 
will be stronger since electrons are always interacting in 
a difference from the case of short range interaction. As 
we will see later the interaction effects will play an im- 
portant role even at low density when the electrons are 
far from each other (i? S> h) and where the interaction 
can lead to the delocalization transition similar to one in 
the 3D Anderson model. It is convenient to study this 
transition by the means of level spacing statistics as it 
was done for 3D one-particle case in ||l5| . 

To analyze the effect of Coulomb interaction between 
two electrons let us consider the 2D Anderson model 
with the diagonal disorder {-W/2 < Ei < W/2), hop- 
ping V, the lattice constant a — 1 and the interaction 
C//|ri — r2|. In these notations Vg — U/{2V^Trns) and 
it is convenient to introduce another dimensionless pa- 
rameter rp = Uli/2^/nV which is equal to r^ value at 
Us = 1/^1. We will consider the case with U ^ V and 
r^ 3> 1 when the average distance between electrons 
i?, = |ri — r2| is much larger than their noninteracting 
localization length: R ~ Xj ^fnl '^ r^ ^ Zi ^ 1. In this 
case the two-body interelectron interaction has a dipole- 
dipole form and is of the order of Udd ^ C/Ari Ar2/i?'^ ^ 



Ulf/R^. Indeed, the first two terms in the expansion of 
Coulomb interaction give only mean field corrections to 
one-particle potential and the nontrivial two-body term 
appears only in the second order in the electron dis- 
placements Ari ^ Ar2 '^ h near their initial positions 
ri 2- The matrix element of this dipole-dipole interac- 
tion between localized noninteracting eigenstates is of 
the order of Us ~ J/X) AriAraV-'^/i?^ ~ U/R^. Here 
-0 ~ exp(— |Ari^2|/^i)/^i are localized one-electron states 
and due to localization the sum runs over Ij sites and 
each term in the sum has a random sign. According to 
the Fermi golden rule these matrix elements give the in- 
teraction induced transition rate Fg ^^ U^P2 ^ Uj^^/V , 
where the density of coupled states in the middle of the 
energy band is still p2 ~ li/V since due to localization 
only jumps on a distance li are allowed. These interac- 
tion induced matrix elements mix two-electron states if 
Ke ~ FeP2 > 1, that corrcsponds to i? < li{Uli/VY^^ 
(a similar estimate for electrons in 3D was given in Ref. 
lb). Since Zi ^ 1 the condition R':^ li is still satisfied. 
For Ke > 1 these transitions lead to a diffusion with the 
rate 
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This diffusion expands in an effective 3D space. Indeed, 
the center of mass of two electrons diffuses in 2D lattice 
plane and in addition the electrons diffusively rotate on 
a ring of radius R and width li. The radius of the ring is 
related to the e-e-energy E ~ U / R which remains con- 
stant. Since R '^ li \t takes a long time to make one 
rotation along the ring. As for the 3D Anderson model 
this diffusion becomes delocalized when the hopping is 
larger than the level spacing between directly coupled 
states, namely: 
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rTirs > 1 
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Formally the situation corresponds to a quasi-two dimen- 



sional case with M, 



ef 
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>> 1 parallel 



planes (number of circles of size li in the ring) so that 
the pair localization length l^ jumps from l^ ^ h for 
Kg < 1 to Zc ~ li exp(7rKe?'/ ) ^ h above the transition 
Ke > 1. The transition is sharp and similar to 3D An- 
derson transition when r^ > r^ 3> 1. If electrons would 
be able to move inside the ring then Mef would be even 



larger (Me/ 
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It is important to stress that the parameter Xe which 
determines the delocalization border and measures the 
effective strength of two-body interaction decreases with 
the increase of r^. Apparently, this looks to be against 
the common lore according to which the larger is r^ 
the stronger is the e-e- interaction. The reason of this 
contradiction with (g) is simply due to the fact that 
Tg compares Egg with Ep computed in the absence of 
disorder. In the presence of not very weak disorder 



{ro = Eee/W <C 1 and r^ :^ 1) the one-electron states 
are localized and form the basis of Coulomb glass p6[ . 
In this Coulomb glass phase the e-e-interaction becomes 
weaker and weaker with the growth of average distance 

— 1/2 

between electrons R '^ Us oc r^ in the natural agree- 
ment with (g). The transition border (||) was obtained 
for excited states. However, it is clear that if the interac- 
tion is not able to delocalize the excited states then the 
low energy states will also remain localized since two- 
electron density p2 drops at low energy. In this sense (g) 
determines the upper border for r^. 

To study the delocalization transition (g) the level 
spacing statistics P{s) is determined numerically for dif- 
ferent system sizes L. To follow the transition from 
the localized phase with the Poisson statistics Pp{s) 
to delocalized one with the Wigner-Dyson statistics 
Pwd{s) it is convenient to use the parameter ij = 
/;°(F(s) - PwD(s))ds//;°(Pp(s) - PwDis))ds, where 
So = 0.4729... is the intersection point of Pp{s) and 
Pwois) 0. In this way yy = 1 corresponds to the Pois- 
sonian case, and r]=0 to Pwd{s)- The dependence of 77 
on the one electron energy e = E/2, counted from the 
ground state is shown in Fig. 1 for different disorder W 
and interaction strength U. Usually ND=4000 realiza- 
tions of disorder are used and in addition the average in 
a small energy interval allows to increase the total statis- 
tics for P{s) and t] from NS = 12000 for low energy 
states up to NS = 10^ at high energies with larger den- 
sity of levels. The matrix diagonalization is done in the 
one-electron eigenbasis truncated at high energies that al- 
lowed to study two electron low energy excitations (with 
energy E) at large system sizes L < 24. The periodic 
boundary conditions are used for one-electron states, the 
Coulomb interaction is taken between electrons in one 
cell of size L and with 8 charge images in nearby 8 cells. 
The Coulomb interaction periodic in one cell gave simi- 
lar results. Only the triplet case was considered but the 
singlet case should give similar results P,p^. 

The results of Fig. la show that at fixed interaction and 
strong disorder W/V — 15 the P{s) statistics approaches 
to the Poisson distribution {rj — 1) at large system size 
L and large rg — U L / {2^/2nV) . This means that all 
states are localized. For smaller disorder the situation 
becomes different (Fig. lb,c). While near the ground 
state still r/ ^ 1 for large L, the tendency is inverted 
above some critical energy Cc where 77 — > 0. All curves 
77(e) for different L are crossed in one point in a way 
similar to the 3D Anderson transition studied in [[l5[ . 
This result can be understood in the following way. At 
strong Coulomb interaction U ^ V the excitation energy 
e is related to the distance between electrons R: e ^ 
U/R (similar relation was used in ||l^ for the Coulomb 
glass). At higher e the distance R becomes smaller, the 
interaction is stronger and for e > Ec the delocalization 
border R ^ U/e ~ liVj^ (||) is crossed and the states 
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FIG. 1. Dependence of rj on the rescaled one-electron en- 
ergy e/B (with B = 4:V) for different W , system size L (a — c) 
and interaction strength U (d). For (a — c): the size is 
L = 6(-^);8(full triangle); 10 (o); 12(D); 16 (full diamond); 
20 (*); 24 (x), so that 2.39 < r^ < 9.57, U/V = 2 and 
W/V = 15(a); 10(6); 7(c). For (d) : W/V = 7,L = 16 and 
U/V = 2 (fuU diamond); l(v);0.4(o); 0.2 (full circle); 0.1 
(full trangle). 

become delocalized. Since the distance R is related 
with the two electron energy E — 2e ^ U /R the spacing 
statistics P(s), which is local in energy and therefore also 
in i?,, is not influenced by states where particles are far 
from each other. In this sense the situation is different 
from the case of short range interaction. According to 
the above arguments ic = ej-i / B should remain con- 
stant when li changes with disorder. The value of Zi can 
be extracted from the average inverse participation ra- 
tio fi = l/X^IV'l^ computed for one-particle states in 
the middle of the band {li ^ ^/^)■ For L = 24 and 
W/V = 10; 7; 5 we have respectively 6 = 11-6; 36.7; 84.2 
that with ec/B « 0.6; 0.28; 0.16 (the case W/V = 5 is 
not shown) gives ic = 3.08 ± 0.01 in a satisfactory agree- 
ment with the above expectations. The variation of i] 
with the interaction U is shown in Fig. Id. According 
to it rj increases with the decrease of U (states become 
more localized) in agreement with the general estimate 
(g). The analysis above allows to understand the depen- 
dence of 77 on e and L. Another reason for the decrease of 
77 at higher e is related to the fact that the two-electron 
density of states p2 grows with energy that allows to mix 
levels more easily. A more detail theory should take this 
fact into account but also to analyze the variation of the 
rate Fg with e. The results in this direction will be pub- 
lished elsewhere |1J] . 

The P{s) statistics for two electrons in 2D near the 
critical point Sc/ B is shown in Fig. 2. Its comparison 
with the critical statistics in 3D Anderson model taken 



from ||l^ (see also |^^) demonstrates that both statistics 
are really very close in agreement with the arguments 
given above. At the critical point the value of r]c is close 
to its value in the Anderson model {ijc = 0.20). The 
small deviations from this value in the case of 2D elec- 
trons (77c ~ 0.25{W/V = 10);0.17(VF/F = 7)) can be 

1/3 
attributed to the fact that the parameter l^' was not suf- 
ficiently large. The investigation of the case with larger li 
requires a significant increase of the system size L > 24. 
Indeed, for L = 24 and W/V = 5 the localization length 
becomes comparable with L {li ^ y^ ~ 9) that gives a 
decrease of rjc « 0.13. 
1.0 



Ts oc Tic since near the critical point (0) Kg ^ 1 and 
rg ~ r^ . It is also interesting to remark that the scaling 
index i^ « 1.5 found in Q is close to the index z/ « 1.5 
near 3D Anderson transition (the fact that in 3D v Ri s 
can be related to the observed symmetry of TV curves) . 
I thank Y.Hancin and A. Hamilton for the stimulating 
discussions of experimental results, D.Braun for the pos- 
sibility to use the data of Ref. |19) and K.Frahni for a 
useful suggestion. 




FIG. 2. Level statistics P{s) for two 2D electrons at the 
critical point: {+) - W/V = 10, L = 12 (0.55 < e/B < 0.65), 
total statistics NS = 4 x 10^ (see Fig. lb); (o) - 
W/V = 7,L = 16 (0.25 < e/B < 0.3), iVS' = 5 x 10^ (see 
Fig. Ic). The full line shows the critical P{s) in 3D Ander- 
son model {W/V = 16.5, L = 14, taken from Ref. [19]); the 
dashed lines give Poisson statistics and Wigner surmise. 

Of course, one cannot expect that the simple model of 
two electrons considered above will explain the variety of 
experimental results obtained by different groups |^-||. 
However, it shows some tendencies which are in agree- 
ment with the experiment. Indeed at large r^ (density 
lower than some critical Uc) the experiments demonstrate 
the transition from metal to insulator. According to Fig. 
4 in B the density at the transition ric oc \/ ^frl drops 
exponentially with the increase/decrease of the mobil- 
ity/disorder [i ex 1/IF^. This qualitatively agrees with 
the estimate (0) according to which near the transition 
logric ~ log(l/rs) ~ — logTi ~ — 1/IF^. However, the 
condition r^ 3> r^ seems to be not well satisfied and ap- 
parently multi-electron effects should be also taken into 
account. Another interesting experimental result (Fig. 
2 in [pj) shows that the conductivity Uc near the criti- 
cal point grows with increase of density nr_ or disorder 
W . This is in a qualitative agreement with the estimate 
(0) according to which CTc ^ D^/V ^ l/l\ ex rj^ ex 
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